Gauge-Away Effect in Cold Gases on Optical Lattices 
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It is shown that a simple modification of the geometry in which Raman lasers are applied to a cold gas in 
an optical lattice results in transforming the emerging effective electromagnetic field into a pure gauge. This 
contrived gauge-away effect can be observed experimentally by measuring the Mott-Insulator-to-Superfluid 
critical point. The underlying mechanism for this phenomenon is the ability to engineer the transfer of the 
transverse component of the gauge potential into its longitudinal one. 
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Ultracold gases on optical lattices provide a powerful in- 
strument for creating quantum devices that can simulate, in 
a controlled manner, a variety of condensed matter systems. 
The emerging possibility of acting on the system with care- 
fully tuned electromagnetic fields invites the design of a new 
generation of experiments (]]]. In particular, it is in princi- 
ple possible to implement effective Abelian JH, [H 0] and non- 
Abelian [5] gauge fields coupled to the Hubbard model, which 
describes the physics of an ultracold gas on an optical lattice. 
For the off-lattice implementations see 10,01; for the rotation- 
based approach see Jg]. This would provide a remarkable re- 
alization of man-made effective gauge symmetry. 

To obtain a clear understanding of effective gauge fields on 
optical lattices, simple experiments with neat signatures must 
be devised to prove the purported gauge phenomenon. Note 
that, in this context, the gauge fields are not real ones, nor are 
the particles they are coupled to electrically charged particles. 
The system is just behaving according to a Hamiltonian that 
simulates an effective, non-dynamic gauge interaction. There- 
fore, a precise characterization of what the equivalent of gauge 
effects are must be put forward. 

A word of caution, however, is in order. In the cold atom 
system the expectation value of any hermitian operator, in par- 
ticular the vector potential, is accessible experimentally. Con- 
versely, in actual gauge systems, be it a quantum Hall sam- 
ple or quantum chromodynamics, observables must not only 
be hermitian but they must also be gauge-invariant. Since no 
real dynamic gauge-symmetry is present in cold atom simu- 
lations, the exciting possibility of performing experiments re- 
vealing information that would not be accessible in the real 
world presents itself. For instance, gauge-dependent effective 
vector fields have been observed in Ref. where a Bose- 
Einstein condensate subjected to an optically generated effec- 
tive magnetic field was studied experimentally. The authors 
reported distinguishing between different zero magnetic field 
configurations. Also, see the viewpoint fioll . 

We here propose a scheme to experimentally simulate a 
very specific gauge property of the interactions that are con- 
trolled on an optical lattice. The idea consists in exploiting 
a simple variation of the geometry of a known experimen- 
tal setup to 'gauge-away' the effective, non-dynamic Abelian 
gauge field. That is, while no change is made on the intensity 
of the lasers acting on the optical lattice, the Abelian elec- 
tromagnetic potential they effectively generate can be trans- 
formed at will from a would-be physical field into a would-be 



pure gauge, i.e. with zero electric and magnetic field. As 
a consequence, the measurement of any gauge-invariant ob- 
servable must reflect the disappearance of the effective gauge 
field. Thus, a gauge-away effect could be observed by means 
of a simple angle rotation of external laser fields. 

Another way of understanding our proposal is as follows. 
The effective gauge fields generated in the optical lattices are 
non-dynamic. The only way to have a gauge transformation 
is to modify the action of the external lasers. In a way, we 
can physically perform gauge transformations. We can also 
modify at will whatever part of a gauge field is longitudinal 
or transverse and check how this affects observables. The 
gauge-away effect corresponds to specifying the experimen- 
tal manipulation that produces a transfer from the transverse 
to the longitudinal components of a gauge potential. 

From a theoretical point of view, our modification of the 
setup of cold gases on optical lattices coupled to external mag- 
netic fields needs a detailed analysis. We shall first present 
the basic elements of our proposal so as to describe an ex- 
perimental signature that would reveal the gauge-away effect. 
Then, we shall analyze some of the subtleties that appear in 
the theoretical description of the system due to the change in 
the configuration of the experiment. 

Let us begin our discussion by recalling that cold bosonic 
atoms propagating on a sufficiently deep optical lattice are ad- 
equately described by the single-band Bose-Hubbard model 
Illl,[l2tl . In recent years, several proposals have been put forth 
whereby the standard setup for generating an optical lattice in 
the lab is modified in such a way that the hopping terms in 
the Hubbard Hamiltonian acquire a position-dependent phase. 
For an infinite two-dimensional lattice the Hamiltonian reads 
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where o\ n n and a m „ create and destroy an atom at a lattice 
site (to, n), respectively, and obey the usual bosonic commu- 
tation relations. The constant J x (J y ) is the site-to-site tun- 
neling energy in the x (y) direction. The parameter U is the 



2 



pair interaction energy at each site, \i m>n is the local chemi- 
cal potential, and N, n n is the local occupation. Most remark- 
ably, bosons hopping on the lattice acquire position-dependent 
phases 6 x (m, n) and 9 v (m, n). 

The previous Hamiltonian closely resembles that of 
charged particles on a lattice interacting with a classical mag- 
netic field perpendicular to it. As in the case of ordinary 
charged particles, from Eq. Q] it is clear that bosons travel- 
ing around a closed loop pick up a phase shift proportional to 
the magnetic flux going through the area bounded by the loop. 
Henceforth, we will identify the phases B % m n , with the com- 
ponents of an effective vector potential A e J (m, n), i = x, y. 

Note that A^f {rn,n) is only an effective non-dynamic de- 
scription of the system. It is not a fundamental field, it does 
not propagate, it is not subject to gauge symmetry. 

Nevertheless, we can make physical modifications of ex- 
ternal lasers that mimic gauge transformations. This permits 
the following symmetry. The modified Hubbard Hamiltonian 
is invariant under a local transformation of A e J (m, n) and of 
the creation and destruction operators given by 

A e J \m,n) — > A e J (m,n) - (A(m + l,n) - A(m,n)), 
Ay f (m, n) — > A y *(m, n) — (A(m, n + 1) — A(m, n)), 
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for any function A(m, n), which is vividly reminiscent of 
a gauge transformation in an ordinary electromagnetic sys- 
tem. Let us stress here that, despite the invariance of the 
Hamiltonian under Eq. [2] in the cold atom setup the phases 
9 l mn = A e /(m,n), and consequently the effective gauge 
field, are physical quantities. In order to perform a gauge 
transformation on the effective vector potential one must alter 
the configuration of the experimental setup, as the local phases 
d l mn ~ A^ (m, n) must be changed. Note that the invariance 
of the Hamiltonian under Eq. |2]maps the solutions of gauge- 
equivalent configurations to each other if only gauge-invariant 
observables are considered. This fact is far from obvious from 
the point of view of the standard theory of degenerate gases. 

Let us now present a way to realize the above modified 
Hubbard model that allows for an experimental verification 
of a gauge-away effect. The basic idea remains to change ap- 
propriately the action of external laser fields so as to transfer 
the transverse part of the generated effective gauge field into a 
longitudinal part. Then, gauge-invariant observables will only 
depend on the former one. 

Of the different ways to create an artificial magnetic field in 
a two-dimensional non-rotating optical lattice 01, let us 
focus our attention on the experimental scheme proposed in 
Ref. J2l . In this setup, atoms in different hyperfine states are 
trapped in different columns of the lattice by adequately po- 
larizing the standing wave of light that forms it. The hopping 
in the direction of these columns, x, is unaltered and as usual 
it is controlled only by the lattice depth in that direction, V x . 
The hopping in the other direction, y, is controlled by extra 
lasers which induce Raman transitions between the two hy- 
perfine states with a complex Rabi frequency. When an atom 



in a given hyperfine state changes to the other state, it is no 
longer at a minimum of the lattice potential and is compelled 
to move to one of the neighboring columns. In order to ensure 
that the tunneling along y is solely controlled by these addi- 
tional lasers, the relation V x << V y must hold. This allows 
us to set the non-Raman induced hopping in the y direction, 
J y , safely to zero in the same way J z must be negligible for a 
two-dimensional optical lattice. 

Let \g) and e) be the two hyperfine levels. The Rabi fre- 
quency at the lattice position x = (m, n/2)a is f2(x) = 
fioe ±l< ' kc ~ ks ^ x , where k 9 and k e are the wave-vectors of the 
Raman lasers, £Iq is a constant and a is the lattice spacing. As 
pointed out in the original work, the tunneling amplitude in the 
y direction is a complex function and is related to the average 
of the Rabi frequency f2(x) with Wannier functions, localized 
and orthogonal superpositions of Bloch waves, at consecutive 
sites. Contrary to the original proposal we will not restrict the 
vector q = k e — k g to lie along the x direction. This modifi- 
cation is what allows the hopping phase in the y direction to 
depend on m as well as on n, 
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where q is the modulus of q and 6 is the angle between q 
and the .T-axis of the lattice. The proportionality constant be- 
tween Ay* (m, n) and (k e — k g ) • x will be computed later on. 
In this setup, the Raman transitions induce no translations in 
the x-axis direction, hence A e J is identically zero for all val- 
ues of 6 and q. Moreover the n-dependence of A e J is a pure 
gauge, i.e. it has no effect on the simulated magnetic field or 
on any other gauge-invariant observable. That is to say, for 
= 7r/2, we have a non-trivial, position-dependent hopping 
phase A y *(n) — ^qa, while the magnetic field is null, 



B z = A ef (m + l,n) - A y f (m 7 n) 
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Consequently, by rotating the Raman lasers to the special an- 
gle 8 = 7r/2 the simulated magnetic field B z is gauged-away 
and the system is in effect simulating a trivial field in a highly 
non-trivial way. 

From a practical point of view, gauge-invariant observables 
can be computed in any gauge. Therefore, we are free to 
shift A y f — > Ay* to a more convenient gauge. To be precise, 
we can work out gauge-invariant observables using the gauge 
transformation A(n) = qasmO / An(n — 1) that removes the 
?i-dependence of A y f . Then, for all gauge-invariant compu- 
tational purposes, the effective Hamiltonian that governs the 
bosons on the optical lattice is 
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where J is calculated by taking the modulus of the average of 
f2(x) with Wannier functions centered at consecutive sites and 
now Ay (m) = ma cos 9, where a 
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FIG. 1: Diagram depicting the geometry of the configuration. The 
Raman lasers, with wave vectors k 9 and k e , are rotated simultane- 
ously with respect to the lattice in such a way that the angle 9 between 
the difference q = k e — k g and the x-axis, and thus the effective 
magnetic flux per plaquette, changes. The pure gauge configuration 
corresponds to 9 = n/2, that is, for q parallel to the y-axis. 



Working with the transformed Hamiltonian, the gauge- 
away effect can be rephrased. In the configuration 9 = ir/2, 
i.e. A e J — 0, the system must be oblivious to the fact that the 
hopping in the y direction carries the original associated phase 
ffjh n . In other words, due to the fact that for this special angle 
the system is gauge-equivalent to the Bose-Hubbard Hamilto- 
nian with real, constant hopping amplitudes it follows that the 
actual system in the lab must behave accordingly. 

To realize the gauge-away effect experimentally it is first 
necessary to identify an observable. A natural candidate 
arises in the context of the Mott insulator (MI) to superfluid 
(SF) transition, which has been studied extensively both the- 
oretically and experimentally in cold atoms on optical lat- 
tices (see Ref. Jgl] and references therein). The experi- 
mental signature of the MI phase is the absence of structure 
(peaks) in the momentum distribution of the gas, G(k) ~ 

u '(k)| 2 Srr' e * fc '^ R_R '' ) ( a R a R')' where u>(k) stands for the 
Wannier function in momentum space. Conversely, in the SF 
phase the momentum distribution shows pronounced peaks at 
the momentum values of the reciprocal lattice. This is due to 
the fact that for the MI phase, | (a R a R> ) | — > 0, while in the SF 
phase (a R a R >) — > constant, for |R — R'| — > oo. 

A crucial observation is that |(<%a R »)| at large separation 
acts as a gauge-invariant order parameter for the MI-SF tran- 
sition. Indeed, (a R a R> ) transforms multiplicatively with the 
phase e l ( A ( R - )-A ( R '' ) - ) and it only vanishes in the SF phase. It 
follows that the MI-SF phase diagram in the J/U — a plain 
can serve as a witness to the gauge-away effect. 

The exact computation of the order parameter for the MI- 
SF transition is not possible due to the difficulty of solving 
the modified Hubbard model exactly for a large lattice. This 
observable can be approximated by a mean field computation, 
which is known to reflect the MI-SF phase transition. To com- 
pute an approximate ground state of the Hamiltonian in Eq. [5] 
and determine the critical point at different values of the mag- 
netic flux, we proceed by decoupling the Hilbert spaces of ad- 
jacent lattice sites in occupation space by the replacement lfl3tl 

a m+l,n a ™,™ * ^m,n a m+l,n ' i l , m+l.n a m,n ~ VVn.nVVre+l.n- 
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FIG. 2: Critical parameter (J/U) c as a function of the magnetic flux 
<E> = a cos 9. The phase diagram was determined using the mean 
field approximation for a 40 x 40 square lattice with open boundary 
conditions for fj, m ,n = 0.5 and with no trap. As the angle of ex- 
ternal lasers is changed, the transfer from transverse to longitudinal 
effective gauge field makes the MI-SF transition to change the lo- 
cation of its critical point. This corresponds to a gauge-away effect 
of the effective gauge field. For a cos 9 — > 0, the hopping phase is 
9 y (m, n) = Ay' (m, n) — an/2, and the MI-SF transition point ap- 
pears at the known value (J/U) c = .043. Although site-dependent 
phases are present, this transition point remains unchanged because 
(m, n) is completely gauged-away, as it only depends on n. 
Only the interval <3> £ [0, 0.5] is considered because the Hamilto- 
nian in Eq. [5] is invariant under the replacement $ — > 1 — $. 



The quantities ?p m ,n, computed according to the relation 
ipm,n = ( a ln+i n)> aie variational parameters over which the 
energy is to be minimized. A random initial set of ip m . n al- 
lows for a diagonalization of the Hamiltonian in the truncated 
occupation basis. The eigenvector with lowest lying eigen- 
value is then used to compute a new set of ip m ,n- The process 
is then repeated recursively and stopped when V>m,n converges 
to a fixed point. 

The numerical convergence of the iterative mean field 
method is severely affected by the addition of the effective 
gauge field. We find many different minimum energy solu- 
tions of i/>m,iv This is consistent with the fact that the Hamil- 
tonian commutes with the momentum in the y direction and 
the ground state is degenerate. In addition, for certain val- 
ues of a cos 9 and for certain initial conditions the process 
does not converge, it oscillates between two different non- 
vanishing values of of i\) m ,n- Nevertheless, this does not occur 
in the MI phase, where ip m ,n = is the only attractor, so it 
should not affect the phase diagram appreciably. The results 
are consistent with those obtained in a perturbation expansion 
in J (3, and with those obtained by a modified mean field 
approach in Ref. lfl5ll . 

The phase diagram we obtain numerically, see Fig. [2] 
shows an approximately linear dependence of (J/U) c on 
a cos 9 for a range of values of cos 9. As the total effective 
magnetic flux per unit plaquette, $ = a cos 9, decreases, the 
transverse component of the effective vector potential is trans- 
ferred to its longitudinal component, until it is purely longitu- 
dinal for cos 9 — 0. When this particular configuration is 
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reached the MI-SF transition takes place at the same value 
of (J/U) c as it does for the un-modified Hubbard model be- 
cause it is a gauge-invariant observable that does not depend 
on the longitudinal part of the vector potential. That is, the 
position-dependent hopping phase introduced by the Raman 
lasers 6 v (m, n) = an/2 has no effect (it is gauged away) be- 
cause it is identified with the longitudinal component of the 
the gauge field, Ay (m,n) = an/2. In summary, the gauge 
away effect we present can then be experimentally observed 
as an increase of the point at which the MI-SF phase transition 
takes places when the angle of Raman lasers is varied. 

Let us now discuss two technical details that have been 
left apart. First, we need to check that the modified Hub- 
bard model is a good representation of the system and, second, 
the parameters of the model can be fixed at our convenience. 
We start by recalling that the matrix element corresponding 
to the Rabi oscillation between the \e) and Jo) states takes the 
form of an average over Wannier functions y2|], J™ m (m, n) = 
| J d 2 xw*(x — x m . n ) fi w(x — x m . n+ i). Using the fact that 
the Wannier functions factorize, w(x) = w x (x)w y (y), it is 

jRam = |n oe iA;'(m,n) J & x ' \ Wx ,( x ')\ 2 cos (2acos9 x') 

x J dy'w;,(y')c^ asi » e y^w y ,(y'-^), (6) 

where Ay — 2ira(cos 9 m + sin#n/2). The difference be- 
tween Ay and the actual gauge field Ay is given by a spa- 
tially constant (but 9 dependent) phase contribution coming 
from the integral in the second line and can be omitted. 

We are ready to address the first technicality. For Eq. [6] 
to be a reliable calculation of the Raman-assisted tunnelling 
it is necessary that be small with respect to the optical 
potential that forms the lattice. Furthermore, the tilting of the 
lattice, a crucial ingredient in this construction, must in turn be 
small compared to the energy gap to the second band and be 
large compared to Hq. A numerical evaluation of J = | J^ am \ 
shows all the above relations can be satisfied simultaneously 
for all 9 while meeting the requirement J = J x . 

The second technicality is perhaps a more relevant issue. 
The experimental setup must remain well calibrated for any 
value of the magnetic flux < <& < 1. From Eq. [6] 
J = J(a, 9), so the ratio J / J x could significantly vary in 
<3>. If this were the case the degree of anisotropy in the model 



would be ^-dependent and could obscure the effects of the 
effective magnetic field. For different choices of the lat- 
tice potentials Vi, i — x,y,z compatible with tight-binding 
approximation and with the requirement that the free hop- 
ping J y and J z be suppressed shows that this problem can 
be avoided by working at fixed a. For instance, in the case 
Vox = 16.Br, Vo y = 25Er, where En is the lattice's recoil 
energy, the maximum variation J in a is more than 25% while 
the maximum variation J in 9 is 4%. Note that the two vari- 
ations go in opposite directions. It follows that it is possi- 
ble, at least in principle, to define the curves <5($) and 6>(<I>) 
such that J remain precisely constant, J(a(<I>), #($)) = J x , 
< $ < 1. 

The following experimental procedure could be used in or- 
der to minimize the deviation of J / J x from unity. Once suit- 
able values of the Vi have been chosen, we start with a pure 
gauge configuration (i.e. 9 = 7r/2). To maximize a, we fix the 
relative angle of the lasers to tt. At this point, we can perform 
diffusion experiments as in lfl2ll for different Rabi frequen- 
cies f^o (whose value is controlled by the frequencies of the 
Raman lasers). By determining the value of Qo for which dif- 
fusion rates in the x and y directions become equal, we can, 
at the same time, check that the pure gauge configuration cor- 
responds to the free hopping model with J (9 = tt/2) = J x 
and prepare the system for the addition of the effective ex- 
ternal magnetic field. Indeed, to turn on B z we only need to 
change the angle 9 from tt/2 without modifying any of the 
other parameters. 

In conclusion, we have shown that the Hubbard Hamilto- 
nian with position-dependent complex hopping phases serves 
as a laboratory to create man-made effective gauge transfor- 
mations. It is then possible to devise a gauge-away experi- 
ment, based on the idea of transferring all the transverse part 
of the effective gauge potential to its longitudinal one, only 
by modifying the angle between the Raman laser fields and 
the optical lattice. Thus, a gauge-invariant observable such as 
the critical value (J/U) c for the MI-SF phase transition point 
changes in a predictable manner as the effective vector poten- 
tial is gauged away. 
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